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1.0 SUMMARY 


A Mach/CAS control law was designed and developed for use on the NASA Terminal 
System Research Vehicle (TSRV) B737 aircraft to conduct research in profile descent proce- 
dures and approach energy management. The system operates by using the elevator to control 
speed; the throttle position being fixed. 

The control law was designed primarily using linear analysis techniques, although a nonlin- 
ear feedforward term was included in the system to improve response to a ramp input. The 
system was modeled on a nonlinear aircraft simulator to confirm the original linear design 
and validate the system design at additional flight conditions. 

The system satisfied all the design requirements for gain and phase margin, transient speed 
target overshoot, speed holding in wind shear, and elevator activity in turbulence. 
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2.0 INTRODUCTION 


This report documents the design and performance analysis of a Mach/CAS control law 
for implementation on the NASA TSRV B737 aircraft. The law is designed to enable NASA to 
perform approach energy management research. For example, on starting descent to an air- 
port, the pilot switches the autopilot from a straight level flight mode to a constant Mach 
descent mode. As the aircraft descends, air density increases and thus CAS increases. At some 
preselected value of CAS, the autopilot switches from Mach hold mode to CAS hold mode and 
continues descending. The Mach, CAS commands, and the switching point may have been 
determined by an energy management computer to satisfy target objectives of distance versus 
altitude, speed, and time of arrival. 

The system operates by using the elevator to control speed. The throttle normally would be set 
at idle thrust, although other settings may be selected to allow for operational considerations 
(e.g., pressurization or anti-icing considerations). 

The report details the linear design results and describes modifications made to the basic 
design to improve performance for ramp input demands and in wind shear conditions. Consid- 
eration is given to a practical implementation yielding transient free switching between Mach 
and CAS mode and also to minimize modification to the existing aircraft. 

The preferred design configuration was implemented on a nonlinear aircraft simulator at the 
Boeing Renton Flight Simulation Center. The nonlinear simulation was used to confirm the 
original design results and validate the design at additional flight conditions. 


3.0 SYMBOLS 


K 


•V2 


a 


VCAS 

VcASg 

®cmd 

^CAScmD 

MachcMD 

^TcmD 

5e 

®c 

F 

Mach 

r 

<l>u 


KUg 


*^Wg 

Ve 

VcLIM 

CJf 

V 

’ wo 

7 

K„K„Kf 


V 


GT 


V 

V 


Actuator gain (scheduled) 

Angle of attack (deg) 

Calibrated airspeed (kn) 

Calibrated airspeed error (kn) 

Co mma nded attitude 
Commanded calibrated airspeed (kn) 
Commanded Mach number 
Commanded true airspeed (kn) 

Commanded elevator deflection 
Conversion factor (CAS to TAS) 

Current measured Mach number 
Damping ratio 

Dryden wind spectrum (longitudinal gusting) 
Dryden wind spectrum (vertical gusting) 
Elevator deflection (deg) 

Equivalent airspeed (kn) 

Feedforward rate limit 
Filtered pitch rate (deg/s) 

Filtered velocity (kn) 

Flightpath angle 
Gains (fixed) 

Groundspeed (kn) 

Initial condition for wind velocity (kn) 
Inertial acceleration of airplane (kn/s) 
Inertial velocity of airplane (kn) 


K, 


X 


Ug 

H 

h 


q,e 


<^u 

a 

Machj^ifjjj 

Vt 


^TcmD 



u 

w 

VmW 
Fhws 
Lu, Lw 

^wind 


ACSL 

ATOPS 

CAS 

RFSL 

TAS 


Integrator gain 
Longitudinal gust (m/s) 

Nominal height (m) 

Perturbed height (m) 

Pitch rate (rad/s) 

Pitch angle (rad) 

RMS longitudinal gust level (m/s) 

RMS vertical gust velocity (m/s) 

Speed of sound 

Switch velocity (CAS to Mach) 

Switch velocity (Mach to CAS) (kn) 

T im a constant for complementary filter (sec) 

True airspeed (kn) 

True airspeed command (kn) 

True airspeed error 
Velocity error gain 

Velocity of airplane in longitudinal axis 
Vertical gust (m/s) 

Wind shear magnitude (kn) 

Wind shear scaling factor (horizontal) (m) 

Wind turbulence characteristic lengths (longitudinal and vertical) 
Wind velocity (kn) 

ACRONYMS 

Advanced continuous simulation language 
Advanced transport operating systems 
Calibrated airspeed 
Renton Flight Simulation Laboratory 
True airspeed 
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4.0 DESIGN AND PERFORMANCE REQUIREMENTS 


T 3 ensure satisfactory performance of the final control law, the following design and 
performance requirements were imposed on the design. The control law must: 

1. Be suitable for use at all flight conditions and throttle settings. 

2. Be implemented with minimum modification to the existing autopilot design (refs. 1 
through 4) commensurate with meeting all stability and performance requirements. 

3. Maintain approximately constant transient response performance at all flight condi- 
tions with a speed control bandwidth equal to or greater than 0.05 r/s. 

4. Provide satisfactory damping (approximately f = 0.7) on dominant roots at all flight 
conditions and demonstrate insensitivity in damping and bandwidth to 10% varia- 
tions in nominal gains. 

5. Limit overshoots in indicated airspeed to within 2 kn for step inputs; maintain indi- 
cated airspeed to within 2 kn for a ramp input of 1 fl/s^. For wind shear of 1 kn/s the 
maximum error in airspeed should not exceed 10 kn and the steady state error should 
be zero. 

6. Maintain gain and phase margins equal to or greater than 6 dB and 60 deg, 
respectively. 

7. Elevator activity should not be greater than comparable existing control systems. 

The Mach/CAS control system can be considered to consist of two control loops (fig. 1): 

1. An outer velocity hold loop using true airspeed (Vj) and acceleration (v) feedback. 

The reference input, either Mach or CAS command, is converted to a true airspeed command 
that drives the outer (velocity hold) loop. As the Mach/CAS conversion is outside the control 
loop, the linear analysis is restricted to considering the attitude loop and the true airspeed 
control loop. 

The reason for converting to true airspeed instead of controlling CAS or Mach directly is that 
this approach requires only one control law. Also, true airspeed is compatible with inertial 
acceleration feedback, which is used to provide damping of the velocity loop. This way, approx- 
imately constant inner loop dynamics are achieved over the entire flight envelope without the 
need for outer loop gain programmers. 

The practical implementation of generating the true airspeed command from either CAS or 
Mach is considered in a later section of this report. 

The design of the control law was accomplished by first considering the inner attitude loop in 
isolation (i.e., d and q feedback only). An attitude control loop was designed that gave suitable 
performance over the whole flight regime. Next, the velocity loop was added and the perfor- 
mance of the overall loop optimized. 
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4.1 INNER AUTOPILOT LOOP 


The inner autopilot loop (fig. 2) was modeled using a gain (KV 2 ) and a first-order lag at 20 
rad/s to represent the elevator actuator. The q feedback was obtained from the baseline control 
system and was prefiltered through a 16 rad/s washout filter and scaled in deg/s. The washout 
filter originally was added to eliminate steady state biases on pitch rate feedback. 

The inner autopilot loop was designed to give approximately constant performance in terms of 
bandwidth and transient response. The three cases selected for examination were considered 
representative of the whole flight envelope. They consisted of: 

1. High speed: CAS = 320 kn, H = 6100m, flaps = 0 deg 

2. Mid speed: CAS = 250 kn, H = 3048m, flaps = 0 deg 

3. Low speed: CAS = 120 kn, H = 500m, flaps = 40 deg 

In the current aircraft design, the forward path gain Kvj is preprogrammed to vary with Vcas 
from 1.0 at 120 kn to 0.275 at 360 kn (fig. 3). It was considered desirable to maintain this gain 
scheduling. 

Root locus techniques were used to examine the low-speed and high-speed performance and 
determine acceptable values for K, and K,. Table 1 shows the dominant complex roots for the 
nominal low-speed gain values K,, = 8 and Kq = 5.16. 


Speed 

w(r/s) 

f 

Low (120 kn) 

9.27 

0.98 

Mid (250 kn) 

16.0 

0.6 

High (320 kn) 

15.9 

0.6 


Table 1. Inner Loop Dominant Complex Roots for Nominal Gains 


The root loci for the low-speed case is shown in Figure 4 for gain variations of K# and Kq 
around their nominal values. It can be seen that the system is well damped for large varia- 
tions in gains K„ Kq ( ± 30%). The corresponding root locus plot at high speed (fig. 5) shows 
lower damping (f = 0.6); however, the system is not sensitive to gain variation of ± 10%. The 
transient response at high and low speed (for a step change in pitch command, dj^j) was 
considered satisfactory (fig. 6). The rise time (time of 63% of final value) was 0.8 sec compared 
with 0.75 sec at low speed. 

A summary of the design results for the inner autopilot loop is shown in Table 2. The table 
includes the design results achieved for the mid-speed case. The minimum gain and phase 
margins were 23 dB and 82 deg respectively for the three cases considered. A typical open loop 
frequency response is shown in Figure 7 (low-speed condition). 
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K, 

m 

Settling 
time (sec) 
to 95% 

(r/s) 

Gain 

margin 

(dB) 

Phase 

margin 

(deg) 



High speed (320 kn) 


5.16 

8 

0.80 

3.05 


4.5 

23.5 

95.2 

Mid speed (250 kn) 


5.16 

8 

0.75 

2.6 


4.8 

23.3 

91.1 

Low speed (120 kn) 

1.0 

5.16 

8 

0.75 

2.25 


3.0 

23.3 

82.0 


Table 2. Inner Autopilot Loop Performance 


4.2 VELOCITY HOLD LOOP 

The basic velocity hold loop (fig. 8) was designed with the same objectives as for the inner 
autopilot loop (i.e., to minimize gain variation and give similar transient performance over the 
whole flight regime). A summary of the results for the three speed cases is given in Table 3 for 
Kj = 8, Kq = 5.16, Kvj = 0.1, K, = 0.35. The dominant complex roots are shown in Thble 4. 



Rise 

time (sec) 
to 63% 

Settling 
time (sec) 
to 95% 

(r/s) 

Gain 

margin 

(dB) 

Phase 

margin 

(deg) 

(3 dB) 

(90 deg) 

High speed (320 kn) 

12.5 

19.5 

0.175 

0.142 

13.8 

65.5 

Mid speed (250 kn) 

12.5 

19.5 

0.18 

0.15 

12.2 

66.3 

Low speed (120 kn) 

15.0 

39.0 

0.09 

0.13 

18.9 

75.9 


Tbble 3. Outer Loop Performance 


Speed 

o)(r/s) 

r 

Low (120 kn) 

0.33 

0.72 

Mid (250 kn) 

0.25 

0.62 

High (320 kn) 

0.23 

0.70 


Table 4. Outer Velocity Hold Loop Dominant Oomplex Roots for Nominal Gains 


Considering the low-speed case the rise time (to 63% of final value) was 15.0 sec (fig. 9). 
Adequate gain and phase margins of 18.9 dB and 75.9 deg, respectively, were obtained (fig. 10) 
and the system closed loop bandwidth (to the 90 deg phase shift point) was 0.16 r/s (fig. 11). 
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The root locus plots show the system sensitivity to gain variations of 50%, and it can be seen 
that the damping of the dominant complex pole (speed control) is sufficiently insensitive to 
gain variations (figs. 12 and 13). 

In the high-speed case, the rise time to 63% had decreased to 12.5 sec (fig. 14). Stability 
margins were again more than adequate at 13.8 dB and 65.5 deg, respectively (fig. 15). The 
root loci of the dominant complex poles are shown in Figure 16 together with the line of 
constant damping for f = 0.7. The diagrams show the insensitivity of the system dynamics to 
gain variation. 

Results at mid speed were very similar to the high-speed case for both the transient and 
frequency response (table 4). 

It was considered desirable to restrict the maximum vertical acceleration to 0.1 g„. It can be 
seen from Figure 17 that a 1-kn step velocity command (high-speed case) gave a maximvun 
vertical acceleration of -0.32 kn/s (-0.017 g„). Hence, a 6-kn change in input would just exceed 
the 0.1 g„ criterion. A limiter on airspeed error was included to overcome this problem. 
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5.0 DESIGN OF NONLINEAR FEEDFORWARD TERM 


The control law described previously, gave good performance to commanded step changes 
in velocity. However, it resulted in poor performance for a ramp input (i.e., constant accelera- 
tion demand). A 0.6 kn/s (1 ft/s) raimp input produced a steady state velocity error of 6 kn (fig. 
18), whereas the design requirement allows 2 kn for this condition. 

This problem could be solved by using one of two methods: 

1. Additional forward path integrator. 

2. Nonlinear feedforward term. 

The second method, the nonlinear term, avoided the possibility of stability problems and was 
considered a simpler option to implement practically. 

The feedforward term consisted of essentially a differentiator followed by a rate limited lag. 
For implementation purposes, the differentiator was modeled by a gain and a high-pass filter 
(wsishout filter) having a small time constant compared with the bandwidth of the velocity 
hold loop. As the dominant loop time constant was 5 sec, a time constant at least ten times 
smaller was required and a value of 0.1 sec was nominally selected. This derivative term can 
be added to the control system error signal and can be interpreted as a shaped command 
input. 

Consider the action of the derivative term alone (fig. 19), that is, one without any rate limiter 
or lag, then for a ramp input, the output from the differentiator is a step. By selection of the 
appropriate gain, the addition of this signal to the loop error signal drives the system to null 
the velocity error due to the ramp input (fig. 20). 

However, for a step input, the output of the differentiator is an impulse, and the addition of 
this signal to the loop error signal causes an undesirable overshoot in the system response 
(fig. 21). 

A solution to the problem is the introduction of a nonlinear element, namely a lag having a 
rate limit (Vcum) 22). This has the effect of severely attenuating the feedforward impulse 
response for a step command in velocity. It can be seen in Figure 23 that lowering the limit 
reduces the overshoot on the step response such that for V^lim s 0.6 kn/s, the response ap- 
proximates the step response for the baseline system without the feedforward term. 

With a ramp input, the output from the differentiator is a step, and the rate limited lag only 
affects the settling time of the lag. For Kp = 10.8 and washout filter time constant of 0.1, the 
feedforward gain is approximately 1. Therefore, for a ramp input of 0.6 kn/s, the limit only 
affects the ramp response for value of Vclim < 0.6 kn/s. This is verified in Figure 24 which 
shows the same ramp response for ^ 0.6 and Vclim = 0.06. 
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The results presented above were carried out at high (Vcas = 320 kn, H = 6100m) and the 
gain Kp was optimized to null the error in TAS (Kp = 10.8). However, at low speed the opti- 
mum response was obtained with the gain Kp = 15. This is due to the fact that for the same 
input (whether step or ramp), the response of the feedback terms and v is not identical over 
the whole flight envelope, whereas, obviously the feedforward term is the same. Rather than 
have gain scheduling, a compromise value of gain was selected (Kp = 13.0) that gave accept- 
able performance over the whole aerodynamic range. Maximum error to ramp input of 0.6 
kn/s was within ±2 kn for both the high- and low-speed cases (figs. 25 and 26). 
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6.0 WIND SHEAR AND TURBULENCE 


The effect of wind shear and turbulence on the control law was investigated using tran- 
sient and frequency response techniques. 

For the control system in Figure 8, a wind shear of 1 kn/s gave rise to 10 kn steady state error 
in airspeed (fig. 27). This system used inertial acceleration to provide the derivative feedback 
of velocity, (jood tracking of wind shear without the standoff error could be obtained by feed- 
ing back the derivatives of airspeed. However, this would result in excessive elevator activity 
in turbulent situations. The addition of a complementary filter in the derivative feedback path 
allows a compromise to be made between the maximum transient airspeed error in a wind 
shear situation and the elevator activity experienced in turbulence. The complementary filter 
consists of: 

(1) Differentiating amd low-pass filtering the airspeed signal. This gives long-term deriv- 
ative feedback of airspeed and eliminates the standoff error. The turbulence compon- 
ent is filtered to reduce elevator activity. 

(2) High-pass (or washout) filtering the inertial acceleration signal and adding the signal 
to the airspeed signal. The two signals complement each other giving a unity transfer 
function. The high-frequency component of the inertial signal is comparatively free of 
noise. 

The implementation of this filter is shown in Figure 28. The maximum airspeed error in a 
wind shear and the time required for the error in airspeed to reduce to zero are a function of 
the time constant Ty. However, the penalty for reducing the effect of wind shear is an increase 
in the elevator activity due to turbulence. The relationships between elevator activity and 
maximum CAS error are shown in Figure 29 for a Dryden turbulence spectrum (ref. 5) with an 
rms level of 0.3 m/s (combined horizontal and vertical). 

As the control law is not to be used for final approach, then the performance criterion in wind 
shear is not severe (i.e., 10 kn max error for 1 kn/s wind shear). The second constraint consid- 
ered was that of elevator activity, which should not be worse than results obtained with the 
existing TSRV aircraft. As the airplane does not have a velocity hold mode that uses the 
elevator, then results were taken from turbulence results for the TSRV glideslope control law. 
The airplane does have this mode but not using elevator. Although the results were not di- 
rectly comparable, they provided a guide to the levels of elevator activity that were 
acceptable. 

From Figure 29, an initial value of xy = 10 sec was taken. For this system, the effect of 
turbulence was determined by modeling the wind using the Dryden spectrum for horizontal 
and vertical and evaluating the rms elevator activity. The results of vertical and longitudinal 
turbulence of 1.5 m/s rms are shown in Ikble 5 below. When compared with the results for the 
glideslope track law, the level of elevator activity was considered acceptable. 
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■ 

Low-speed elevator 

Activity (deg rms) 


Velocity hold law 

Glideslope track law 

Vertical gusts (1 .5 m/s rms) 

1.4 

1.56 

Longitudinal gusts (1 .5 m/s rms) 

0.5 

0.835 


Table 5. The rms Elevator Activity Due to Turbulence 


The effect of the complementary filter (tv = 10 sec) on the frequency response of the elevator in 
longitudinal turbulence can be seen in Figure 30. Above approximately 0.3 r/s, the comple- 
mentary filter raises the frequency response by approximately 7 dB; however, this was consid- 
ered acceptable. At low speed (fig. 31), the increase in the frequency response caused by the 
addition of the filter was not significant above 0.1 r/s. However, below 0.03 r/s the filter 
reduced the sensitivity to gust, i.e., improved performance to wind shear. 

In the case of vertical turbulence, the effect on the elevator activity of including the comple- 
mentary filter was not significant (fig. 32) for both the high and low speed. 

The effect of the complementary filter on the response of airspeed error for longitudinal turbu- 
lence can be seen in Figure 33. The filter has the effect of reducing the sensitivity to wind 
shear (i.e., low-frequency sensitivity) without significantly increasing the sensitivity to 
turbulence. 

The sensitivity of the airspeed error to vertical turbulence was not significantly affected by 
the inclusion of the filter in both the high- and low-speed cases (fig. 34). 


7.0 IMPLEMENTATION OF MACH/CAS HOLD 


7.1 CONVERSION TO TRUE AIRSPEED 


Consideration has been given to the method of converting from CAS and Mach to true air- 
speed, which is used in the control law. Figure 1 shows a simplified block diagram of the 
control system using a true airspeed command (Vtcmd) airspeed (V<p) feedback to form 

the error signal 


= Vi 


( 1 ) 


It is possible to compute from the pilot selected Vcasc^d ® conversion formula. 

However, the conversion formula is complex and requires accurate measurement of several 
parameters. The risk of incurring a significant error in is therefore substantial. Such 

an error would cause an equivalent standoff between the selected Vcascmd f^® actual 
VcAs- This approach is therefore unacceptable. However, this problem can be eliminated by 
implementing the control loop in terms of the error in Vcas- 


^CASg ^CAS(;j^(ij * ^CAS 

and conversion of this error signal into the equivalent true airspeed error signal 


(2) 


Vte - F.Vc«^ (3) 

In this case, a less accurate conversion factor (derived from speed conversion tables) will 
suffice. 


F = 1 (4) 

(1 - 0.394 X 10- •*/i) 

where h is in m. As is controlled to zero, so Vcas^ must g® to zero, and an exact measure- 
ment of h is not critical to accurately controlling Vcas- 


In the Mach hold control loop, the accuracy problem can be avoided. Restating Equation 1 
above 


Vte = Vtomd - Vt 

then can be obtained from = 

where a = speed of sound and MachcMO = 


MachcMD • a 

commanded Mach number. 


therefore 


„ _ MachcMn • Vt 

t'cmd Mach 

(5) 

where Mach = Current Mach number. 


V = Machg^p y y 

Mach ^ ^ 

(6) 

V = /Machp^p j 

V Mach J ^ 

(7) 


It can be seen that errors in measuring Vt will not cause steady state errors in Mach hold. As 
Mach tends to MachcMo» so tends to zero. The implementation of the equations discussed 
above is shown in Figure 35. 
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7.2 CONTROL LAW INITIAL CONDITIONS 


During the first section of descent, the control system is designed to hold the aircraft at 
constant Mach. When CAS equals a preselected CAS, then the system automatically switches 
to CAS hold and continues the descent by maintaining constant CAS. On ascent, the reverse is 
true. The control law initially is designed to maintain constant CAS, and at a preselected 
Mach, the system switches to holding constant Mach. 

In order to m inimiz e transients on switching between the two modes of operation, it is neces- 
S8uy to set 

Vcascmd » (8) 

when in the Mach hold mode (in descent), and switch to CAS hold when 

VcAs “ ^CASswrrcH • 

Similarly, when the airplane is climbing and in the CAS mode, the Mach command should be 
set to the current value, i.e.: 

MachcMD = Mach (10) 

and the system switched to Mach hold when 

Mach * Machg^^rrrcH (11) 

An additional important consideration is that the forward loop integrator should be held at its 
current value (clamped) if the elevator limits during flight. The integrator should be “un- 
damped” when the elevator comes of the limit or the input to the integrator changes sign. 
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8.0 VALIDATION OF RESULTS 


The equations describing the Mach/CAS control law were implemented on a nonlinear 
aircraft simulator. This was considered desirable in order to: 

1. Validate the design results obtained using linear analysis design techniques and 
results obtained using a simplified nonlinear model with linear aerodynamics (ACSL 
Simulation (ref. 6)). 

2. Examine the noise free performance at other flight conditions not covered by the 
linear analysis. 

3. Demonstrate the ability to provide transient free switching from Mach hold to CAS 
hold and vice versa. 

4. Confirm the performance of the system in wind shear and turbulent conditions. 

8.1 NOISE FREE RESULTS 

Initially, the free airplane responses obtained using the ACSL program were compared with 
the results obtained with nonlinear aircraft simulation. 

The ACSL simulation used derivatives obtained at two flight conditions namely: 

1. Low speed: Vc^s = 120 kn, H = 500m, 7 = 0, and flaps = 40 deg 

2. High speed: Vc^s = 320 kn, H = 6100m, 7 = 0, and flaps = 0 deg 

Simulation runs were carried out at these conditions with an initial offset of 1 deg from 
trimmed condition. Figiire 36a shows plots of 6 and a for the linear simulation at the low- 
speed conditions and Figure 36b shows the nonlinear simulator results for 6 and a. 

It can be seen that the ACSL simulation gave a phugoid oscillation of 37 sec compared with 42 
sec with the nonlinear simulator. The short period response of a crossed the trimmed steady 
state value at 1.6 sec and 4.0 sec (for the linear system) compared with 2.0 sec and 4.5 sec for 
the nonlinear aircraft simulation. 

At high-speed, the linear system took 106 sec (fig. 37a) for the phugoid oscillation compared 
with 82 sec (fig. 37b) for the nonlinear aircraft simulation. The short period crossing points (a 
response) occurred at 0.64 and 1.8 for the ACSL system and 0.7 and 2.2 for the nonlinear 
aircraft simulation model. 

A detailed comparison of the B-737 aerodynamics on the nonlinear aircraft simulation was not 
possible within the time scale available. However, the discrepancy between the ACSL and 
nonlinear aircraft simulation was acceptable considering the difference in the levels of sophis- 
tication between the two models. 

To confirm the validity of the overall control law, the time response to step change in true 
airspeed was examined. Figure 38a shows that for the linear model at high speed, it took 20 
sec to achieve 95% of the steady state value compared with 18 sec for the nonlinear aircraft 
simulation result. At low speed, the 95% level was achieved after 36.8 sec using the ACSL 
simulation (fig. 38b) and 35 sec with the nonlinear aircraft simulation. This small discrepancy 
in the models was considered acceptable and could be explained by differences in the aerody- 
namics and by the additional complexity and nonlinear elements of the nonlinear aircraft 
simulation. 
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small discrepancy in the models was considered acceptable and could be explained by differ- 
ences in the aerodynamics and by the additional complexity and nonlinear elements of the 
nonlinear aircraft simulation. 

A disadvantage of the linear simulation was its restriction to a limited number of flight 
conditions. The addition of extra flight conditions would require extensive modification to the 
program. However, the nonlinear aircraft simulation allowed the performance of the system to 
be determined at other flight conditions, namely; 

1. CAS = 350 kn, H = 6100m, flaps 0 deg 

2. CAS = 250 kn, H = 6100m, flaps 0 deg 

3. CAS = 200 kn, H = 3048m, flaps 0 deg 

4. CAS = 160 kn, H = 500m, flaps 15 deg 

The response to a step change in true airspeed at these flight conditions is shown in Figure 39. 
It can be seen that the performance is within the design requirements for all conditions. 

In addition to the step response, the effect of a ramp input (0.3 m/s^) was examined. In order to 
achieve a similar response at both high and low speed compared with the ACSL results, a 
feedforward gain (Kp) of 6.5 was selected (figs. 40 and 41) compared with linear design require- 
ment of 13.0 (figs. 25 and 26). 

A primary requirement of the system is that it should provide transient free switching be- 
tween the Mach and CAS hold mode. This is demonstrated in Figure 42, which shows an 
aircraft descent from approximately 6700m, at a rate of -7.6 m/s. The Mach hold mode has 
been selected so that initially the system captvires and holds Mach = 0.7. When holding Mach 
number and descending, CAS increased. The simulation stopped when CAS was 324.7 kn, 
although the measured value tended toward the requirement of 325 kn. Maximum transient 
error at switchover was less than 1 kn. 

Transient free switching between CAS and Mach numbers was also demonstrated. For exam- 
ple, Figure 43 shows the effect of Mach and CAS for an aircraft climbing from 6100m at an 
initial nominal climb rate of 7.6 ra/s. Initially, the control system is in the CAS hold mode, 
maintaining 300 kn while Mach number increases. At Mach = 0.68, the Mach hold mode 
automatically is switched in and the system holds Mach with a transient at switchover of less 
than Mach 0.003. 

8.2 PERFORMANCE IN WIND SHEAR AND TURBULENCE 

Figure 44 shows the results of the ACSL simulation for a horizontal wind shear of 1 kn/s 
plotted in terms of Mach response. Peak error in terms of airspeed was 7 kn and the error was 
reduced to 10% of its peak value by 38 sec. Figure 44 also shows the results of the nonlinear 
simulation, which gives approximately the same peak error and settling time as the ACSL 
simulation. 

To determine the effect of wind shear at low speed, a wind shear profile was selected. The wind 
velocity is given by the equation; 

Vmw = FHWS • Vmwic 

where FHWS is a scaling factor dependent on height (fig. 45) and V^vvic is an initial condition 
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for wind velocity set to 1 kn. The velocity error (VcASg) shown in Figure 46 for an aircraft 
climbing at a nominal 7.6 m/s. It can be seen that the maximum error in CAS was 0.32 kn. 


In addition to wind shear, the effect of turbulence was examined at high and low speed. A 
Dryden wind spectrum (ref. 5) was selected, which is defined as 

^ 2 a„^Lu 

<t>. (y) = ^ 


((1+(L„7)^) 

for longitudinal gusting, and 

^ (-y) = g.r'Lw(l + 3(L^ 

’• T (1 +(LW7)^)^ 

where Lu, Lw are the characteristic lengths 

a„, a, are the rms gust velocity 

and 7 - w/v. 


Typical values of scale length and rms level were selected, i.e., Lu, Lw = 183m; <t„, = 1.5m 

rms. A t 5 rpical plot of gust is shown in Figure 47; the rms level for the 100 sec sample was 
1.4m. The effect of turbulence on the elevator is shown in Figure 48. 


A summary of the results is shown in Table 6. For a medium level of turbulence (1.5m rms) the 
elevator response was considered acceptable. 



Longitudinal 

Vertical 

Inputs rms (m/s) 

1.38* 

1.38* 

Elevator rms response (deg) 

0.26 

0.7 

Velocity rms response (kn) 

2.83 

0.47 


• Calculated for 100 sec time history. 

Table 6. Summary of Turbulence Results 



9.0 CONCLUSIONS 


1. A Mach/CAS hold control system was designed using linear analysis techniques for 
flight conditions. 

• High speed: CAS * 320 kn, H = 6100m, 7 = 0, flaps = 0 deg 

• Mid speed: CAS = 250 kn, H = 3048m, 7 = 0, flaps = 0 deg 

• Low speed: CAS = 120 kn, H = 500m, 7 = 0, flaps = 40 deg 

2. The basic design was modified— 

• lb improve performance in wind shear by the addition of a complementary filter in 
the feedback path. 

• To improve the response to a ramp input by the inclusion of a nonlinear feedforward 
term. 

2. The basic design was modified— 

• To improve performance in wind shear by the addition of a complementary filter in 
the feedback path. 

• To improve the response to a ramp input by the inclusion of a nonlinear feedforward 
term. 

3. For the design cases considered, the system satisfied the design criteria. 

• The minimum gain and phase margins were 12 dB and 65 deg. 

• The rise time (to 63% of final value) for a step input was 13.5 sec ± 11%. 

• The overshoot to a step input was less than 2 kn. 

• The peak error in airspeed, due to a 1 kn/s wind shear, was less than 10 kn. 

The error reduced to zero in the steady state. 

• Elevator activity was less than 1.5 deg rms for vertical gusts of 1.5 m/s rms. 

• For a ramp input of 1 kn/s^, the error in airspeed was less than 2 kn. 

4. Consideration was given to practical implementation in order to provide transient free 
switching between the Mach and CAS hold modes. The system has been designed to be 
implemented with minimum modification to the existing design. 

5. The control system was implemented on the nonlinear aircraft simulator, which enabled 
simulation runs to be made at several flight conditions not covered in the linear design, 
namely: 

• CAS = 160 kn, H = 500m, flaps = 15 deg 

• CAS = 200 kn, H = 3048m, flaps = 0 deg 

• CAS = 250 kn, H = 6100m, flaps = 0 deg 

• CAS = 350 kn, H = 6100m, flaps = 0 deg 

6. The nonlinear aircraft simulation runs confirmed that the system satisfied the design 
requirements for noise free, wind shear, and turbulent conditions. 
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Figure 5. Root Locus of Dominant Compiex Poles Inner Loop, High-Speed Case 
(6, 6 Feedback Only) 
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Figure 7. Low-Speed, Open Loop Frequency Response 
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Figure 10. Velocity Hold Loop Frequency Response (Low Speed), Open Loop 
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Figure 13. Root Locus for Dominant Complex Poles (Outer Loop, Low-Speed Case) 
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Figure 21. Derivative Feedforward lerm: Step Response 
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Figure 22. Mach/CAS Control Law With Feedforward Term With Rate Limited Lag 
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Control system as shown in Figure 22 
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Figure 24. Ramp Response for ^nation in Rate Limit 
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High speed 



Figure 26. Final High-Speed Response 
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Figure 28. Mach/CAS Control Law With Complementary Filter in Velocity Feedback Loop 



CAS error due to wind shear of 1 kn/s 
Elevator activity due to turbulence of 1 ft/s2 
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Dryden turbulence spectrum 



Figure 29. Elevator Activity and CAS Error Versus Complementary Filter Time 
Constant for Wind Shear of 1 kn/s 


48 



Filter 


No filter 


Frequency, rad/s 

Figure 30. Frequency Response of Elevator (deg) (With and Without Complementary 
Filter) (High Speed) ^ 
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Figure 31. Frequency Response of Elevator 6 q. (deg) (With and Without Complementary 
Filter) (Low Speed) ^ 
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Figure 32. Elevator Frequertcy Response in Vertical Turbulence 
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Figure 34. Airspeed Error Frequency Response in Verticai, Turbulence 








(a) ACSL simulation 



(b) Nonlinear simulation 
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Figure 36. Transient Response of Free Aircraft (1-deg initiai Offset on a and 6) (Low Speed) 
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(b) Nonlinear simulation 



Figure 37. Transient Response of Free Aircraft (1-deg initiai Offset on a and 6) 
(High Speed) 
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Figure 38. Step Response of Velocity Hold Loop 











Figure 41. Ramp Response High Speed (RFSL Simulation) 
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Figure 42. Mach Hold to CAS Hold Aircraft Descent 
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Figure 43. Ascent: CAS Hold to Mach Hold 
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Figure 46. Effect of Wind Shear on Velocity (Low Speed) 
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